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Abstract. Observations of magnetically active close binaries with orbital periods of a few days reveal the existence of starspots 
at preferred longitudes (with respect to the direction of the companion star). We numerically investigate the non-linear dynamics 
and evolution of magnetic flux tubes in the convection zone of a fast-rotating component of a close binary system and explore 
whether the tidal effects are able to generate non-uniformities in the surface distribution of erupting flux tubes. Assuming a 
synchronised system with a rotation period of two days and consisting of two solar-type components, both the tidal force and 
the deviation of the stellar structure from spherical shape are considered in lowest-order perturbation theory. The magnetic field 
is initially stored in the form of toroidal magnetic flux rings within the stably stratified overshoot region beneath the convection 
zone. Once the field has grown sufficiently strong, instabilities initiate the formation of rising flux loops, which rise through 
the convection zone and emerge at the stellar surface. We find that although the magnitude of tidal effects is rather small, they 
nevertheless lead to the formation of clusters of flux tube eruptions at preferred longitudes on opposite sides of the star, which 
result from the cumulative and resonant character of the action of tidal effects on rising flux tubes. The longitude distribution 
of the clusters depends on the initial parameters of flux tubes in the overshoot region like magnetic field strength and latitude, 
implying that there is no globally unique preferred longitude along a fixed direction. 
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1. Introduction ing non-linear numerical simulations we investigate the influ- 

^, . , , . ,. , . , r ■ , • cnce of tidal effects on the dynamics and evolution of mag- 

Observational results indicate that in close, rast-rotating bi- . „ , ■ ■ j , ■ ■, , 

^, , . T^^r ... netic nux tubes inside the convection zone until they erupt 

naries like RS CVn and BY Dra systems magnetic activity , .,1. n c r. 1 • j . 

„ at the stellar surface. Particular attention is paid to the non- 

3 IS usually much more vigorous than in the case of the Sun .„ . . . ^, , . „ j- . 1. r 

-^1 — "^l 1 uniformities in the resulting surface distributions of erupting 

" " (see. e.g. JStrassmeier et al.ll993L and references therein). Data a . , j .u • j j • 1 j-.- -ru- 

" ' ■ • flux tubes and their dependence on initial conditions. This 

sets covering time periods much longer than the rotation pe- , , ^ . ^ ,. • . 

^ ^ , .... .. . work complements our previous studies of the equilibrium and 

nod of a system reveal non-uniform longitudinal spot distri- ,. , ... ■ ^ ■ r, , ■ , 

^ : , ■ s linear stability properties of magnetic flux tubes in the over- 

butions (with respect to the direction of the companion star), . , • . 1 1. j r 

^ r /' shoot region at the low er boundary of the convection zone 

which indicate the existence of long-lasting spot conglomer- 1,, , \ o o i,-- — n-irvrvn '^nnnfi. nc^r, 

, ^ , (IHolzwarth & Schussleii2 000. 200_3|, hereafter HSI). 
ations at preferred longitudes where starspots persist longer 

or flux eruption is more frequent than at other surface re- 
gions ( Henrv et alJ 119 95: Jetsu 1996). The preferred longi- 
tudes are frequently about 1 80° apart on opposite sides of Sect. [^contains a brief summary of the magnetic flux tube 
the stellar hemisphere (e.g., IBerdvugina & Tuominenlll998t model and the analytical approximation of the close binary 
IBerdvugina et al.lll998l) and their dependence on rotation pe- model. In Sect.|3]we explain the accomplishment of the numer- 
riod, i.e., essentially on the distance between the two compo- ical simulations as weU as the implementation of initial condi- 
nents, suggests that they are related to t he proximity of the tions. Sect. |3 shows the resulting latitudinal, longitudinal and 
companion star taeckert & Ordw ay'l995V surface distributions of erupting flux tubes in the case of a bi- 
In analogy to the 'solar paradigm' for magnetic activity, nary systems with an orbital period of two days as well as a 
it is assumed that starspots at the stellar surface are gener- study concerning the dependence of non-uniformities on the 
ated by erupting magnetic flux tubes which originate from the orbital period. Sect. |5] contains a discussion of our results and 
lower part of the outer convection zone of the star. By us- Sect. |6]gives our conclusions. 
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2. Model assumptions 

2.1. The flux tube model 

We study the spot distribution on active binary stars by ap- 
plying the 'solar paradigm' for magnetic activity, assuming 
that starspots are due to erupting magnetic flux tubes (e.g., 
ISchiissler et alJ I199(tI) . In this model, the magnetic field is 
stored in the form of toroidal flux tubes inside the subadia- 
batic overshoot-region below the convection zone at the inter- 
face to the radiative core. Once a critical magnetic field strength 
is exceeded, perturbations initiate the onset of an undulatory 
(Parker-type) instability and the growth of flux loops, which 
rise through the convection zone and lead to the formation of 
bipolar active regions and starspots upon emergence at the stel- 
lar surface. 

The numerical code we have used for the non-linear simula- 
tionsis explicated in Moreno-Insertis ( 1986) and Caligari et al. 
ill99.4 . All calculations are carri ed out in the framework of 
the 'thin flux tube approximation' (ISpruiil 98lS describing the 
evolution of an isolated magnetic flux tube embedded in a field- 
free environment. The surrounding convection zone is treated 
as an ideal plasma with vanishing viscosity and infinite con- 
ductivity, implying the concept of frozen-in magnetic fields. 
The dynamical interaction between the flux tube and its envi- 
ronment is taken into account by an hydrodynamic drag force 
acting upon a tube with circular cross-section which moves rel- 
ative to its environment. The radius of the cross-section is de- 
termined according to the amount of magnetic flux and the lat- 
eral balance between the external gas pressure and the internal 
total (gas and magnetic) pressure. Since the radiative time scale 
in the convection zone typically exceeds the rising time of flux 
tubes, we assume the adiabatic evolution of flux tubes. 



2.2. The binary model 

We consider a detached binary system with two solar-type 
components, which move on circular orbits with an orbital pe- 
riod of two days. For the active star we assume synchronised, 
solid-body rotation, with the spin axis orientated perpendicu- 
lar to the orbital plane. Its internal structure is described by a 
tidally deformed solar model, whereas the companion star is 
considered as a point mass. The model includes both the tidal 
force and the deviation of the stellar structure from spherical 
shape, which are sufficiently small to be treated in lowest-order 
perturbation theory (see HS/, Sect. 2.2). 

Using spherical coordinates (radius r, azimuth 0, latitude 
A; (p = Q is the direction toward the companion star), the dom- 
inating term of, e.g., the azimuthal component of the effective 
gravitation, g-gff , can be written as 
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where g-^ is the gravitation in the unperturbed case, q is the 
mass ratio between the two components, e = (r/a) the ex- 
pansion parameter of the approximation and a the distance be- 
tween the two components. For a weakly deformed star it can 
be assumed that the value of any stellar quantity, /, (pressure, 
density,. . . ) is still only a function of the effective potential, 



*Peff - "P* + *Ptide, wherc "P* is the gravitational potential of the 
star and ^Ptide the potential of the tidal perturbation [see Eq. (2) 
in HS/]. Following Eq. (5) in HS/, the (Eulerian) perturbation 
of the stellar structure. A/, can be written as 
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where the undisturbed value, /*, and its local scale height, ///, 
are determined from an unperturbed stellar model. The leading 
orders of the tidal effects given in Eqs. Q and (|2j, respectively, 
both show a 7r-periodicity in the azimuthal direction (p. 

Using Kepler's third law, the order of magnitude of tidal 
effects becomes 
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while owing to the ratio r/Hf in Eq. the perturbation of 
the stellar structure also depends on the scale height of the con- 
sidered quantity. It is enhanced if / changes in a relative small 
stellar layer, like the superadiabaticity, (5 = V - Vad, inside the 
thin overshoot region below the convection zone. For the sys- 
tem considered here (T = 2d, = Mq, q - \ and r < Rq), 
Eq. Q yields e^q ~ 10"^ (with a ~ 8 Rq). 

3. Simulations 

The linear stability analysis in HS/ describes the response of 
a toroidal equilibrium flux tube to small perturbations inside 
the overshoot region and its evolution prior to its penetration 
to the convection zone above. In the evolution of rapidly ris- 
ing flux loops in the superadiabatical region non-linear ef- 
fects of the tube dynamics become important, which have to 
be followed by numerical simulations. According to the re- 
sults in HS/, we consider flux tubes with initial magnetic field 
strengths between 7 ■ lO'^G < B, < 2 ■ 10^ G. The initial 
flux rings are located in the middle of the overshoot region 
at latitudes between 0° < /Is < 80° parallel to the equato- 
rial plane and Uable to an undulatory (Parker-type) instability 
(e.g.. lFerriz-Mas (feSchiissler 199 5). The simulations start with 
a perturbation of a flux ring and cover up to 32 yrs of its evolu- 
tion in the convection zone. In order to sample the influence of 
tidal effects on rising flux tubes, we carry out a large number 
of simulations in which each tube is subject to a perturbation 
localised around a different longitude, (p^- This is realised by 
an in-phase superposition of sinusoidal displacements, sin nup, 
with azimuthal wave numbers m = 1 ... 5 and amplitudes of 
one percent of the local pressure scale height (Fig.Q. For simu- 
lations pertaining to a given (/?s, /is) configuration, this specific 
perturbation is successively shifted in longitude; the interval 
< 0s < 180° is sampled with 20 simulations, each separated 
by A/Is - 9°. Owing to the azimuthal 7r-periodicity of the un- 
derlying problem, the resulting pattern of flux tube emergence 
is then extended to the interval 180° < 0s < 360°. The regularly 
distributed positions (/Is, 0s) of the localised perturbations are 
shown as points in the polar projection of the overshoot region 
in the left panel of Fig.|2l The evolution of rising flux tubes is 
followed close to their eruption at the stellar surface, resulting 
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Fig. 1. Schematic illustration of the initial 'localised perturba- 
tion' of a flux ring (dashed) resulting from the in-phase super- 
position of sinusoidal wave modes with m = 1 ... 5. The per- 
turbation is largest around the longitude 0s- The central arrow 
indicates the direction toward the companion star. 
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Fig. 2. Polar projections showing the overshoot region with 
the uniformly distributed positions (/Is, 0s) of initial perturba- 
tions {left) and a surface distribution with the positions (/Ig, 0e) 
of erupting flux tubes resulting from their non-linear evolu- 
tion in the convection zone (right). In this example with Z?s = 
1.2 ■ 10'' G preferred longitudes are discernable at intermediate 
latitudes /Ig ~ 45° near the quadrature longitudes; a number 
of flux tubes do not emerge due to their stable initial configu- 
ration. The star rotates in counter-clockwise direction and the 
companion star is located in the direction (p - 0. 



in spot distributions (Ag, 0e) like the example shown in the right 
panel of Fig.|2l A complete set of simulations represents a non- 
linear mapping from the bottom of the convection zone to the 
top, assuming uniform distributions for the perturbations over 
longitude and for the flux rings over latitude in the overshoot 
region. The simulations stop at about 0.98 Rq, where the 'thin 
flux tube approximation' ceases to be valid. Since the rise of 
the tube's crest through the uppermost parts of the convection 
zone is very fast and nearly radial, there is hardly a difference 
between its final position in the simulation and the position of 
eruption at the stellar surface. 

A toroidal isentropic flux tube embedded in the tidally dis- 
torted stratification of a binary star actually represents a non- 
equilibrium configuration, since it exhibits a density contrast 
along the tube axis which entails a small deviation from me- 
chanical equilibrium. Following the linear stabiUty analysis in 
HS/, the dynamical tube evolution due to this azimuthal density 
variation proceeds on very long time scales of several thousand 
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Fig. 3. Eruption times fg in years, where the position (B^,As) 
in the diagram refers to the initial flux tube configuration in the 
overshoot region. According to the linear stability analysis in 
HS/, grey shaded regions mark domains of Parker-type insta- 
bilities with dominating azimuthal wave numbers m - 1 (dark) 
and m = 2 (light). Configurations in white regions are linearly 
stable, but yield eruptions due to the non-linear flow instability. 

days, which are not relevant for the eruption of magnetic flux at 
the stellar surface.. Since the perturbation of the flux ring due 
to the density variation along the tube is much smaller than the 
perturbation arising from the superimposed displacement (de- 
scribed above) and furthermore identical for all simulations at a 
given latitude, its eff'ect, if any, on the longitudinal distribution 
of erupting flux tubes is expected to be neglectable. 

4. Results 

4.1. Eruption times 

The eruption times, fe, defined as the elapsed time between the 
initial perturbation of a flux tube and its emergence at the steflar 
surface, are found to be of the order of several months to years, 
in a few cases (near the instability boundaries) up to decades. 
Their dependence on the initial (B^, A^) configurations (Fig.|5} 
is roughly proportional to the respective growth times of the 
linear instability in the overshoot region, which are, accord- 
ing to the analysis in HS/, dominated by Parker-type instabili- 
ties with azimuthal wave numbers m - 1 and 2. On the other 
hand, flux tubes in linearly stable initial configurations, e.g., at 
As - 65 . . . 75° or in the stability island at intermediate latitudes 
and high field strengths, are found to become unstable by an- 
other mechanism, which we ascribe to the interaction of a flux 
tubes with the environment due to its hydrodynamic drag once 
the relative tangential flow velocity exceeds the phase velocity 
of a perturbation; a detailed description of this flow instability 
will be given in a subsequent paper. This second-order effect is 
not included in the framework of the linear stability analysis in 
HS/, but nevertheless appears in the course of non-linear sim- 
ulations and leads to the formation of rising flux tubes on time 
scales comparable to those of the Parker-type instability. At 
high magnetic field strengths this instability mechanism con- 
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Fig. 4. Poleward deflection, AA = A^-A^, of erupting flux tubes 
as a function of their initial latitude A^ in the overshoot region. 
Symbols mark diff'erent initial field strengths, B^. 



tributes a substantial fraction to the total number of erupting 
flux tubes. In contrast, simulations with initial configurations 
in the instabiUty islands at low field strengths and /Is ~ 20° 
result in eruption times of several hundred years and thus are 
probably irrelevant. 

The simulated rise of flux tubes through the convection 
zone does not always proceed in a continuous manner. Some 
flux tubes become unstable and form rising loops in the over- 
shoot region, but then pass through a dynamical stage of quasi- 
stability where their radial rise temporarily comes to a halt be- 
fore they resume after some several hundred days their ascent 
through the convection zone toward the surface. 

4.2. Latitudinal distributions 

Whereas sunspots exclusively appear in an equatorial belt be- 
tween about +35° latitude, fast-rotating stars frequently show 
spots also at intermediate, high, and even polar latitudes. In 
the framework of the magnetic flux tube model, this behaviour 
is explained by the influence of the Coriolis force, which acts 
on the internal gas flow along the tube axis and causes a pole- 
ward deflection of the tube ( Schiissler & Solanki 199 2). Figure 
l^shows the latitude excursion, AA = A^-As, i.e., the difiference 
between the latitude of emergence, A^, and the starting latitude 
in the overshoot region, for flux tubes with diff'erent initial mag- 
netic field strengths, B^. The deflection is largest for flux tubes 
starting at low latitudes. It decreases considerably with increas- 
ing As and eventually drops to zero at high latitudes. Flux tubes 
starting above A^ > 60° essentially do not show any deflection. 
The poleward deflection depends only marginally on the ini- 
tial magnetic field strength. At moderate latitudes, AA becomes 
somewhat smaller for increasing B^, because the larger buoy- 
ancy force leads to a more radial rise through the convection 
zone. 




0.0 0.1 0.2 0.3 0.4 0.5 



Fig. 5. Distributions of eruption longitudes, 0e, as a function 
of the initial perturbation longitude, 0,. Crosses mark the re- 
sults of simulations, solid lines the (interpolated) non-linear 
functions 0e(0s) in the binary problem, and dotted lines the 
linear functions in the corresponding single star problem. The 
grey histograms at the left side show the resulting distribu- 
tions of erupting flux tubes binned to 9°-wide longitudinal in- 
tervals (abscissa gives normalised values). The distributions 
show (from top to bottom) a rather flat distribution, a clearly 
preferred longitude, and a favoured longitude range for emer- 
gence, respectively. 

4.3. Longitudinal distributions 

Figure|5]shows examples of longitudinal distributions of erupt- 
ing flux tubes in a binary star Owing to the wave-like char- 
acter of the instability, the longitude of emergence, diff'ers 
from the formation longitude of a rising loop in the overshoot 
region. In a single star, an azimuthal shift of the localised per- 
turbation by an ofl'set A^s yields the same off'set in the longi- 
tude of emergence. Under the action of tidal forces in a de- 
formed binary star, however, the dynamics and evolution of a 
flux tube depend on the value of <ps- The longitudes of emer- 
gence thus exhibit a deviation from the uniform distribution 
for a single star, with the amplitude of the 7r-periodic deviation 
depending on the (B^, A^) values of the initial flux ring. This re- 
sults in a non-uniform longitudinal distribution of erupting flux 
tubes (represented in Fig.|5]by the histograms at the left side of 
each panel), which show clear indications for preferred longi- 
tude ranges of flux eruption: beside cases with only marginal, 
if any, longitudinal non-uniformity, there are cases of highly 
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Fig. 6. Longitudinal distributions of erupting flux tubes with 
high initial latitudes, = 65... 80°. The shading indicates 
the amount of clustering inside 9° -wide longitudinal intervals 
normalised to the total number of flux tubes for each parameter 
pair {Bs, As); the numbers give the maximum values for each 
case. Because of the 7r-periodicity in longitude, only the range 
< (pe/iln) < .5 is shown. White regions indicate stable, non- 
erupting flux tubes. 
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peaked distributions as well as distributions where a broad lon- 
gitudinal interval is preferred. 

4.3.1 . High latitudes (^s - 65 . . . 80°) 

Figure|6lshows the longitudinal patterns of emergence for high 
initial latitudes, based on the corresponding histogram for each 
(Bs, As) configuration. The shading indicates the amount of lon- 
gitudinal clustering of erupting flux tubes. For A^ ~ 80° we 
have Parker-type instabilities with the azimuthal wave number 
m - 1 and almost flat distributions in longitude. For latitudes 
between /Is ~ 65 ... 75° and moderate field strengths consider- 
able asymmetries appear in the form of highly peaked distribu- 
tions or broad preferred intervals. In some cases, up to 45% of 
ah flux tubes of a (Bs, '^s) configuration emerge within a lon- 
gitude interval of only 9° width. Apart from a decreasing ten- 
dency for clustering with increasing field strength, there is no 
concise systematic dependency of the longitude distribution of 
these features on the initial (B^, A^) parameters. 

4.3.2. Intermediate latitudes {A, = 30... 60°) 

The evolution of flux tubes starting at intermediate latitudes are 
significantly affected by tidal effects. Figure |7] shows longitu- 
dinal distributions of erupting flux tubes starting between A^ = 
30 . . . 60° in the overshoot region. In the domain of Parker-type 
instabilities with the azimuthal wave number m = 2 (left to 
the vertical black lines), the concentration and orientation of 
clusters of emerging flux tubes depends systematically on the 



Fig. 7. Longitudinal distributions of erupting flux tubes with 
intermediate initial latitudes /Is = 30 . . . 60°. The vertical black 
lines divide the domains of linear instabilities according to HS/: 
the upper left domain is dominated by Parker-type instabilities 
with m - 2, whereas the region to the right shows the flow in- 
stability. White regions indicate stable, non-erupting flux tubes. 

initial values (Bs, /Is)- The most prominent clusters occur for 
As ~ 35 . . .40° and Bs ~ 8 ... 14 ■ 10"* G. An important as- 
pect is the orientation of preferred longitudes with respect to 
the direction of the companion star a.i (j) — 0. Figure shows 
that clusters basically show up at all orientations depending on 
the initial configuration. Since there is no unique orientation 
for all configurations we use the term preferred longitude in 
relation to the clustering of erupting tubes for an individual 
{Bs,As) configuration. For increasing magnetic field strength 
near the instability threshold, the clusters are shifted toward 
smaller longitudes (retrograde with respect to the rotation of 
the system), whereas for high Bs this trend reverses and the 
clusters shift toward larger longitudes (prograde). The slope of 
this longitudinal shift as a function of field strength depends on 
the starting latitude of the tube. At small As, the slope is rather 
large and an increase ofBs by about 5- 10^ G can result in an 
azimuthal shift of the cluster by 90°. The slope decreases con- 
siderably with increasing starting latitude and approaches zero 
at As ~ 55°. Here the orientation of clusters is apparently inde- 
pendent of the magnetic field strength, thus representing a kind 
of fixed preferred longitudes. However, in this case the clus- 
tering is rather weak. The shift rate seems to be related to the 



6 



Holzwarth & Schiissler: Flux tubes in close binaries II. 



amount of clustering, with more concentrated clusters showing 
larger shift rates. 

4.3.3. Low latitudes {A,^0... 25°) 

Longitudinal distributions of erupting flux tubes starting at 
equatorial latitudes (A^ < 10°) and liable to the Parker-type 
instability with wave number m = 1 do not show significant 
non-uniformities. For = 15 . . .25°, flux tubes become un- 
stable above Bs - 1-4 ■ 10^ G with the dominating wave mode 
m -2. This domain of initial configurations is characterised by 
a particular evolutionary behaviour: flux tubes with the same 
(Bs, /Is) do not exhibit a systematic continuous dependence of 
the emergence longitude on (f)^ like in Fig. |5j but an irregular 
distribution. This aspect of the simulation results is discussed 
in more detail in Sect.|5l 

4.4. Surface distributions 

Figure |8] shows polar projections of the surface distributions 
of erupting flux tubes with different initial field strengths, Z?s, 
combining both longitudinal and latitudinal information with 
the grid of simulations (A0s, A/l,) = (9°, 5°) shown in the left 
panel of Fig. |2l Owing to the poleward deflection and the az- 
imuthal non-uniformities induced by tidal effects, the surface 
distributions deviate considerably form the uniform distribu- 
tion of initial flux tubes at the bottom of the convection zone 
and typically show features on opposite sides of the stellar 
hemisphere resulting from the dominating ;7r-periodicity of tidal 
effects. While some parts of the surface are avoided by erupt- 
ing flux tubes, like the equatorial belt below about 25°, they 
gather in other regions to form clusters of flux eruption. For 
low magnetic field strengths, only a few flux tubes starting from 
intermediate latitudes reach the surface, whereas the majority 
of tubes at other latitudes does not leave the overshoot region 
owing to their stable initial configuration. For increasing field 
strengths more flux tubes become unstable and contribute to 
the overall surface pattern erupting particularly at moderate and 
high latitudes. The flux eruption at high latitudes is mainly due 
to the flow instability, which results in clusters, if any, cover- 
ing broad longitudinal intervals. Flux tubes erupting at mod- 
erate latitudes are due to Parker-type instabilities with wave 
number m = 1 starting at low initial latitudes. These config- 
urations hardly show longitudinal asymmetries. However, flux 
tubes starting at somewhat higher latitudes with m = 2 expe- 
rience a smaller poleward deflection and emerge at moderate 
latitudes, too, showing considerably non-uniform longitudinal 
distributions. At intermediate latitudes and high field strengths 
the extended and less pronounced clusters are due to the flow 
instability. The systematic dependence of preferred longitudes, 
which is discernible in the distributions shown in Fig.0 is less 
evident in the surface patterns, because, owing to the variable 
poleward deflection, flux tubes from different initial latitudes 
emerge in the same surface area. Since the orientation of clus- 
ters in distributions with the same field strength but dififerent 
latitudes A^ are usually not the same, their superposition in the 
surface pattern sometimes leads to broader emergence patterns. 
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Fig. 8. Surface distributions of erupting flux tubes for different 
initial field strengths, Bs. The star rotates in counter-clockwise 
direction, cp - Q is the direction toward the companion star. 
The shading indicates the number of flux tubes erupting in a 
(A0e, A/le) - (9°, 5°)-wide surface area. The values in brackets 
give the relative number of flux tubes which do not erupt at the 
surface due to their stable equilibrium in the overshoot region. 
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Fig. 9. Longitudinal distributions of erupting flux tube for sys- 
tems with orbital periods T - \,2 and 5 d (from left to right). 
Shown are distributions for flux tubes starting at intermediate 
latitudes. 

However, particularly at low field strengths, the existence of 
considerable clusters at preferred longitudes remain, while for 
higher field strengths significant non-uniformities are still dis- 
cernible in the overall surface patterns. 

4.5. Dependence on system period 

We have investigated the influence of the magnitude of the tidal 
effects by considering systems with orbital periods T - Id and 
5d, which correspond to binary separations a = 5.3 Rq and 
15.5 Rq, respectively. Figure|9]shows a comparison of longitu- 
dinal distributions of eruption for flux tubes with intermediate 
initial latitudes and field strengths Bs = 1 ... 1.5 ■ 10^ G. In 
the case T - 5d, most longitudinal distributions show hardly 
any significant clusters of eruption. Merely for /Is = 60° there 
are some non-uniformities, which indicate that tidal effects can 
still be effective for some initial tube configurations. The case 
T - 2d, the reference case discussed in detail above, shows the 
existence of considerable emergence clusters and the system- 
atic dependence of their orientation on the initial tube configu- 
ration. For the shorter orbital period T = 1 d the degree of clus- 
tering increases. However, for some initial tube configurations 
there is a transition from highly peaked distributions toward 
more extended preferred intervals. In contrast to highly peaked 
clusters these patterns of eruption are rather characterised by 
an avoidance of certain longitudinal intervals. 

5. Discussion 

Our results show that the presence of the companion star in a 
close binary system can considerably alter the evolution and 
surface distribution of erupting magnetic flux tubes. The actual 
influence of tidal effects on a rising flux tube depends on its ini- 
tial magnetic field strength, B^, and latitude, A^. This includes 
the stability properties of the initial flux ring in the overshoot 
region, i.e., the instability mechanism and the azimuthal wave 
number, m, of the dominating unstable eigenmode, which gov- 



erns the character of the subsequent non-linear evolution. In 
our range of parameters, only single-loop tubes (with m - 1) 
and double-loop tubes (with m - 2) appear. 

The properties of erupting flux tubes can roughly be clas- 
sified according to their linear instability mechanism as deter- 
mined in HS/. In the domain of Parker-type instabilities, erupt- 
ing single-loop tubes starting at either low or very high latitudes 
tend to form uniform longitudinal patterns of eruption, whereas 
double-loop tubes result in clusters of eruption concentrated in 
narrow longitudinal intervals. This specific behaviour is pre- 
sumably due to the congruency between the double-loop ge- 
ometry of the rising tube and the 7r-periodicity of the tidal ef- 
fects, which enables a resonant interaction. Although there is 
no globally unique orientation for the clusters of eruption, the 
longitudinal shift of the clusters (usually prograde with increas- 
ing Bs) can become very small and even vanish leading to fixed 
preferred longitudes at a certain latitude. Tubes erupting due to 
the flow instability show a strong susceptibility to tidal effects 
leading to extended preferred intervals. In contrast to the do- 
main of the Parker-type instability with m - 2, the orientations 
of the preferred intervals in this case lack a systematic depen- 
dence on initial tube parameters. Particularly at high latitudes 
and high field strengths the flow instability is responsible for 
the majority of erupting flux tubes. The dependence of the evo- 
lution of flux tubes on the instability mechanism (with similar 
initial configurations otherwise) indicates that the formation of 
non-uniformities at the surface is not only a matter of the tidal 
forces acting during the non-linear rise of the flux tubes through 
the convection zone, but also affected by stability properties in 
the overshoot region. While the Parker-type instability is ba- 
sically a buoyancy-driven mechanism and thus susceptible to 
perturbations of gravity and stellar structure, the flow instabil- 
ity is more related to the tube parameters, namely its diameter 
and the relative velocity perpendicular to the environment. 

For some (Z?s,/ls) combinations the set of erupting flux 
tubes does not yield a continuous function 0e(0.s) like the ones 
shown in Fig. [S] Instead, they yield either piecewise steady 
functions, which frequently alternate between two branches 
or irregular distributions, which show no obvious relation to 
the original shift of perturbation longitudes, 0s- These partic- 
ular configurations with irregular distributions are limited to 
Bs ^ 1-5 ■ 10^ G, /Is - 15 . . . 35°, a parameter domain mainly 
showing a Parker-type instability with m = 2, which roughly 
coincides with the region of retarded flux eruption described 
in Sect. 14.11 Because of the dominating unstable wave mode 
with m - 2, two loops of similar size form at opposite sides of 
the initially toroidal flux tube in the overshoot region and enter 
the convection zone, whereby the wave-like character of the 
instability causes a longitudinal propagation. Non-linearities 
lead to slightly different evolutions of the two loops, so that 
the early stage of the development is characterised by consec- 
utive changes of the dominance between both loops until one 
of them prevails and eventually emerges. In this stage, the time 
delays of rising flux loops described above entail longitudinal 
offsets owing to their wave-like propagation. Flux tubes start- 
ing with different values A0s may suffer different time delays 
before they resume their rise and emerge at the surface. The 
originally continuous phase shift A(ps imprinted on a set of flux 
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tubes is thus annihilated and the resulting longitudinal distribu- 
tion irregular. Since for any (Z?s, As) configuration only twenty 
simulations (per semi-hemisphere) are accomplished, the irreg- 
ular distributions can however exhibit non-uniformities due to 
the small statistical sample. For a large number of simulations 
we would expect a rather uniform distribution. 

The trajectories of rising flux tubes exhibit the poleward 
deflection typical for fast rotators, which is caused by the 
dominance of the C oriolis force over the buoyancy force 
jSchiissler & SolankHflQQZ) . The latitudinal deflection leads 
to an overlap of patterns of erupting flux tubes with differ- 
ent initial latitudes, which results in a superposition of clus- 
ters at different longitudes thus broadening the range of pre- 
ferred longitudes in the surface distributions. The deflection is 
largest for flux tubes starting at equatorial and low latitudes 
and decreases for intermediate and high initial latitudes. An 
explanation for the observed high-latitude and polar spots in 
the framework of the present erupting flux tube model thus 
requires large amounts of magnetic flux at high latitudes in 
the overshoot region, either generated by a suitable dynamo 
mechanism or transported by large-scale motions inside the 
convection zone. Alternatively, flux erupting at lower latitudes 
may be transported tow ard the pole by meridional motions 
( ISchriiver&Titlel200lh . 

The investigation of binary systems with different orbital 
periods, i.e., distance between the components, shows that the 
non-uniformities become more pronounced with stronger tidal 
effects. The orbital period of about T ~ 5d for the onset of clus- 
tering is in agreement with the results from lHeckert & Ordwavl 
( 1995). In general, the non-uniformities increase for smaller or- 
bital periods, but this trend is accompanied by a transition from 
clusters localised in narrow longitudinal intervals to broader 
intervals, up to major fractions of the stellar circumference. 
The number of cases with irregular longitudinal distributions 
increases for smaller orbital periods. This is in accordance with 
our interpretation that, with a closer proximity of the compan- 
ion star, the resonant interaction between the tube and the tidal 
effects becomes more efficient and eventually ends up in an 
analogy to the 'resonant catastrophe' of a driven oscillator 

The action of tidal effects on the dynamics of magnetic flux 
tubes is basically a cumulative effect, which proceeds and adds 
up during the entire rise from the overshoot region up to the sur- 
face. During this time, an individual flux tube is subject to fur- 
ther interactions with, e.g., turbulent convective motions, dif- 
ferential rotation, meridional flows, and with other concentra- 
tions of magnetic flux. Although these additional interactions 
may modify the evolution and alter the final location of erup- 
tion to some extent, we do not expect them to severely diminish 
the characteristic «on-axisymmetric influence of the compan- 
ion star since they are presumably axisymmetric on average. 
However, it should be kept in mind that the applied thin flux 
tube model represents a rough simplification of the magnetic 
field structure in convective stellar envelopes. 

The distributions of flux eruption shown above allow sta- 
tistical predictions, which apply to large numbers of erupt- 
ing flux tubes with similar initial configurations. Our simu- 
lations are based on flux tubes evenly distributed in latitude 
as well as on perturbations evenly distributed in longitude. 



A changing situation in the overshoot region, e.g., a depen- 
dence of the amount of magnetic flux on latitude due to a 
particular dynamo mechanism or a dependence of tube pertur- 
bations on longitude due to a dependence of convective mo- 
tions on tidal effects, could, in principle, be taken into ac- 
count by folding the specific initial condition with the present 
non-linear mapping of rising flux tubes from the bottom of 
the convection zone to the top. Although the frequently ob- 
served occurrence of spot clusters at opposite sides of stel- 
lar hemispheres can consistently be explained by the domi- 
nance of the TT-periodic contributions of tidal effects, our time- 
independent patterns of emergence do not immediately al- 
low a description of the frequently obser ved azi muthal mi gra- 
tion of spot clusters in time J strassmeie ^|l994t jRodono e t alJ 
1995, 2000; Berdvugina & Tuominen 1998). However, a non- 
stationary dynamo mechanism with latitudinal propagation of 
a dynamo wave would lead to a time- and latitude-dependent 
field strength and magnetic flux, so that the systematic depen- 
dence of the cluster longitudes on field strength shown in Fig.0 
could transform into a longitudinal migration of the spot clus- 
ters in time. 

We considered a binary system with two solar-type main 
sequence components to investigate the basic effects of tidal 
interaction on the evolution of erupting flux tubes. Although 
the results shown here are consistent with some observations, 
there are observed spot properties which are not explained by 
our results, like spots at ver y low latitudes or fixed spo t clusters 
around the substellar point jLanza et alJ200k.200Z:.Olah et alJ 
1^02). RS CVn systems, however, typically consist of at least 
one (sub-)giant component with a much deeper convective en- 
velope than the Sun, leading to the assumption that the forma- 
tion and orientation of spot clusters are not exclusively depen- 
dent on the binary separation alone, but also in a more com- 
plicated way on the type and structure of the star. The deep 
convection zone of a giant, for example, will probably modify 
the properties of patterns of eruption owing to the cumulative 
character of the interaction between the rising flux tube and the 
companion star. 

6. Conclusion 

Erupting flux tubes in binary stars are considerably affected by 
the tidal influence of the companion star. Although the magni- 
tude of tidal effects - the tidal force and the deviation of the 
stellar structure from spherical symmetry - are rather small for 
the systems considered here, they are capable to alter the sur- 
face distributions of erupting flux tubes. This results in clusters 
of emerging flux tubes at preferred longitudes and in the forma- 
tion of spot clusters on opposite sides of the active component. 
We consider two main reasons for this behaviour: 

1. Contrary to single-loop tubes, which are not strongly af- 
fected, the geometry of double-loop tubes (with azimuthal 
wave number m - 2) allows of an efficient resonant inter- 
action with the congruent 7r-periodic structure of the domi- 
nating tidal effects. 

2. Since the rise of magnetic flux tubes from the bottom of the 
convection zone until eruption at the surface takes about 
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several months or years, the cumulation of the tidal effects 
leads to a strong influence on the resulting emergence pat- 
tern. 

The distributions of erupting flux tubes shown here are 
based on uniformly distributed initial conditions (distribu- 
tions of flux tubes and their localised perturbation inside the 
overshoot-region). More realistic spot distributions require the 
application of a dynamo model, which determines the spatial 
and temporal distribution of magnetic flux at the bottom of the 
convection zone. Furthermore, calculated surface distributions 
for one system cannot easily be used for other systems with 
diff'erent stellar or orbital parameters since changes of these 
quantities alter the tube evolution and patterns of emergence; 
individual studies have to be carried out for each system. 
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